Semiclassical Theory of Chaotic Conductors 
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We calculate the Landauer conductance through chaotic ballistic devices in the semiclassical limit, 
to all orders in the inverse number of scattering channels without and with a magnetic field. Families 
of pairs of entrance-to-exit trajectories contribute, similarly to the pairs of periodic orbits making 
up the small-time expansion of the spectral form factor of chaotic dynamics. As a clue to the exact 
result we find that close self-encounters slightly hinder the escape of trajectories into leads. 
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Electron transport through conductors is usually dom- 
inated by scattering from disordered impurities. The 
pertinent theory involves averages over random configu- 
rations of scatterers and ultimately amounts to random- 
matrix theory. Typical consequences of disorder on trans- 
port, like weak localization, actually survive for "clean" 
systems whose linear dimension is smaller than the mean 
free path, provided the ballistic motion is chaotic due 
to irregularly shaped boundaries of the "cavity" serving 
as the conductor [l|, |2| . The equivalence of chaos in the 
ballistic limit and impurity scattering represents a theo- 
retical challenge similar to that of proving the so-called 
Bohigas-Giannoni-Schmit conjecture about universality 
of fluctuations in energy spectra of chaotic dynamics. 

Gutzwiller's semiclassical periodic-orbit theory has re- 
cently been brought to success in explaining spectral fluc- 
tuations of individual chaotic dynamics as faithful to the 
ensemble averages provided by random-matrix theory Q . 
Analogous attempts have been undertaken for ballistic 
transport 0, 0], employing semiclassical techniques in 
Landauer's scattering approach to conductance The 
two lowest-order terms in a l/A^-expansion (with N the 
dimension of the scattering matrix) of the conductance 
were thus found in agreement with RMT. Here we take 
up that line of thought, with two crucial modifications to 
be explained below, which yield the conductance to all 
orders in l/N. Both for systems with and without time 
reversal invariance as well as for the crossover between 
these two cases caused by a magnetic field we confirm 
the equivalence to impurity scattering as predicted by 
random-matrix theory. Each order of the l/N expan- 
sion turns out one-to-one with the small-time expansion 
of the universal spectral form factor given in |3j]. Our 
confirmation of the validity of the RMT conductance for 
individual chaotic cavities with arbitrary N is particu- 
larly relevant for (electronic or microwave) experiments 
involving a number of channels of order unity [(j . 

We imagine a two dimensional cavity connected to the 
outside by two leads respectively supporting N± and N2 
channels. The entrance-to-exit transition amplitudes for 
(electron) waves then form a unitary N x N S'-matrix 
where N = Ni + A^ . Sub-blocks t, t' of respective size 
N1XN2, JV2X Ni comprise the amplitudes for transitions 
between different leads while the remaining sub-blocks 



r, r' refer to "reflections" , i. e. transitions with entrance 
and exit at the same lead. We are concerned with the 
conductance between different leads, given by Landauer 
(3 as G = Tr(tH), in units of the quantum conductance 
Go = e 2 /nh. Random-matrix theory yields an averaged 
conductance Gp, the average being over the Wigner- 
Dyson circular ensembles of unitary matrices from the 
orthogonal {(3 = 1) or unitary (/3 = 2) class @,H|], 
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We first deal with the orthogonal and unitary classes. 
Further below, we shall allow for a magnetic field caus- 
ing a smooth transition between these classes. Random- 
matrix theory provides an integral expression [9] G(£) 
interpolating between G\ and G2, as a parameter £ pro- 
portional to the squared magnetic field is varied. 

We shall derive a semiclassical expression for the con- 
ductance of a single chaotic cavity. The van Vleck prop- 
agator for the amplitudes t brings classical entrance-to- 
exit trajectories into play Q, U t nm oc ^ A 1 e lS i/ h , 
with S* 7 the action and A 1 a stability amplitude of 7; the 
sum includes all trajectories connecting the leads with 
discrete incident angles, one for each channel. The Lan- 
dauer conductance becomes the double sum Q, 



Tr(th) 




a i(^-S v )/B 



(2) 



somewhat similar to Gutzwiller's Q double sum over pe- 
riodic orbits for the spectral form factor; here Tjj is the 
Heisenberg time, related to the volume of the energy shell 
as Tjj = il/2irh. The brackets (•) denote an energy av- 
erage over a small interval near the Fermi energy. 

In the limit fi — > only trajectory pairs with the action 
difference AS* = Sy — Sy of the order h may survive the 
averaging implied in (0). The simplest such family is 
the diagonal one with 7 = 7'. To evaluate the diagonal 
approximation Tr(i t i) d i ag = Tff 1 (Z) 7 = 7 ' IAI 2 )> Richter 
and Sieber invoked a new sum rule which we write as 
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FIG. 1: Pairs of trajectories 7,7' contributing to the conductance: a) a Richter/Sieber pair differing in a 2-encounter, b) a 
pair of trajectories differing in one 2-encounter and one 3-encounter. In a), 7 and 7' differ in the sense of traversal of one loop 
(y = 1), and the relative orientations of stretches in the 2-encounter lead to fi = 0. For b) we have v = 3; y, = — 1 for the 
2-encounter and /1 = I for the 3-encounter 



Like the somewhat analogous sum rule given by Han- 
nay and Ozorio de Almeida 10] for periodic orbits, the 
Richter-Sieber sum reflects the ergodicity of long trajec- 
tories; in addition, it accounts for the escape of a trajec- 
tory hitting a lead after a dwell time T by the exponential 

N rp 

"survival probability" e T « . In spite of its quantum 
guise, which is convenient for our purposes, the escape 
rate N/Th is, of course, a classical quantity; Planck's 
constant cancels since the number of channels is propor- 
tional to the ratio of the width of the leads and the de 
Broglie wavelength. Clearly, the diagonal approximation 
gives the first term of the 1/N expansion of the random- 
matrix conductance and even the full conductance 
for the unitary symmetry class. 

Higher-order terms are due to families of pairs (7,7') 
of trajectories. The partner trajectories in a pair nearly 
coincide, at least up to time reversal, in "loops" whose 
cumulative duration is of the order Th ; these loops are 
separated by close self-encounters of each trajectory; the 
encounter stretches have durations of the order of the 
Ehrenfest time Tg cx lnfi, <C Th oal/h, and connect the 
loops differently for 7 and 7'. The various families of 
pairs (7, 7') can be imagined as arising from the pairs of 
periodic orbits contributing to the analogous expansion 
of the spectral form factor [3 by cutting a loop into two 
and pulling the two open ends to different leads. The 
cutting described actually turns the families of pairs of 
periodic orbits giving the term r™ in the r expansion 
of the spectral form factor into the pairs of trajectories 
contributing to the n-th order in the 1/N expansion of 
the conductance. The dwell time T of 7 and 7' is then the 
cumulative duration of all loops and encounter stretches. 

A clue to the full expansion in search is the survival 

N rp 

probability. That quantity is a little larger than e t h 
for all trajectory pairs except the "trivial" ones, 7' = 7. 
To see that, we focus on, say, an ^-encounter where I 



stretches of 7 (and of 7') are mutually close at least up 
to time reversal and separated from the leads by loops; 
each of the I intra-encounter stretches has the duration 
t C nc- Survival of 7 inside the cavity during the first of 
the I stretches implies survival during the I — 1 other 
stretches. The whole encounter therefore does not con- 
tribute lt cnc but only t enc to the effective exposure time 
in the exponential survival probability. The cumulative 
exposure time is the sum of all loop durations, augmented 
by one instead of I contributions t cnc for each Z-encounter, 

T exp = X/loops ^loop + JJce ^onc — T — ^2 a {la ~ -O^cno 

the index a labeling the encounters. Due to Te <C Tjj 

N rp 

the survival probability e t h exp differs only little from 

N rp 

the naively expected value e T « , but the difference will 
presently be seen to be decisively important. 

We illustrate our procedure for the first correction to 
the diagonal approximation which involves partner tra- 
jectories differing in a single 2-encounter. Fig. la depicts 
the two encounter stretches and three loops with dura- 
tions ti, i = 1,2, 3; the dwell time is T = t\ + t% + t% + 
2t cnc . The two encounter stretches must be nearly mutu- 
ally time reversed since otherwise the partner trajectory 
would decompose such that one of the loops would con- 
nect the two leads while the other would close in on itself 
to form a periodic orbit; this is why the first off-diagonal 
correction requires time reversal invariance. 

We may search for 2-encounters drawing a Poincare 
section V transverse to 7. If the section hits a 2- 
encounter, 7 must pierce through V in two points which 
are almost mutually time reversed. We denote the stable 
and unstable shifts of the second (time reversed) piercing 
relative to the first by s, u. For the quantitative defini- 
tion of 2-encounters we require \u\ < c,\s\ < c with c 
a constant independent of Ti (whose exact value is ir- 
relevant). The coordinates u, s determine the encounter 
duration t cnc = (1/A) lnc 2 /|su| and the action difference 
AS = 5 7 — Sy = us; here A is the Lyapunov constant. 

The ergodic probability of the second piercing oc- 
curring with given s, u and duration ti of the loop 
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separating the encounter stretches is ds du dt 2 /Sl. We 
can write the number of trajectory pairs generated by 
2-encounters with the given action difference AS* as 
dAS Jj s i <c \ u \ <c ds du wt{u, s)6(AS — us). The weight 

function wt(u, s) is obtained from the ergodic return den- 
sity l/fi by (i) integrating over all durations t 2 of the loop 
separating the encounter stretches; (ii) integrating over 
all positions of the Poincare section which is equivalent 
to integrating over the duration of the first loop t\ ; (iii) 
dividing by the encounter duration t cnc (the probability 
to cross an encounter is proportional to its duration; in 
order to obtain the number of encounters from the cross- 
ing probability the encounter duration has to be divided 
out). Requiring all loop times to be positive we get 



W T (u,s) 



dt\dt^ 



e(t 1 )G(t 2 )e(T-t 1 -t 2 -2t enc ) 



(T - 2t cnc (u,s)) 
2t cnc (u, s)Q 



(4) 



The first off-diagonal term in the conductance thus reads 



Each encounter requires its own Poincare section which 
latter contains I close-by piercings by 7; the separations 
of subsequent piercings with respect to the first ones can 
be expressed in terms of coordinates Sj,Uj,j — 1, . . .1 — 1 
obeying |sj|,|«j| < c Q. The ergodic probability den- 
sity for the L — V returns is l/Sl L ~ v . To obtain the 
weight function analogous to Q for the number of en- 
counter sets with a given v, the ergodic density has to 
be integrated over all loop times but one, and divided 
by the product of all encounter durations. In addition 
the result has to be multiplied by the number Af(v) of 
topological versions ("structures") of trajectory pairs dif- 
fering by (i) the order in which encounters are visited in 
7; (ii) the relative direction of stretches in each encounter 
and (iii) the intra-encounter reconnections leading from 
7 to 7'. The action difference in the pair can be written 
as AS = J2 a j s ajUaj- When applying the sum rule we 

use the exposure time T exp = J2f=i *i+E a =i *eno as ex ~ 
plained above. Integrating over the remaining loop time 
th+\ (instead of the dwell time) and summing over all v 
we obtain the complete expansion of the conductance 



G« = (/ dsdue isu / R J2 

\J \s\<c,\u\<c 



I A, 



-w t {s,u) ) . (5) 



The sum over 7 is now calculated using the rule @ with 
the crucial modification pointed out above, i.e. the sur- 
vival probability e ^5 (* 1 +* 2 + t3 +*« ic ). The integrations 
over ti,t2 (in wt) and over the dwell time T can be re- 
placed by integrations over all three loop durations. The 
lower limit of a loop duration is zero which expresses a 
second crucial ingredient of our approach: the encounter 
stretches may neither overlap nor stick out of the cavity 
openings [Tl| : the upper limit is infinity. The result 
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(6) 
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obtains a factor Th/N from each loop-time integration. 
To calculate the integral over s, u we expanded in powers 

— — t N 

of the encounter duration, e t h onc = 1 — ^-i enc + . . .; 
the contribution of the zero-order term is annihilated by 
averaging; terms of second and higher order vanish in 
the limit ?i — > |3| ; in the remaining first-order term t enc 
cancels against the denominator whereupon the integral 
becomes trivial, ± J^ c dsdue ius ^ h -> 2Trh/fl = 1/T H . 

Higher orders in 1/N are due to pairs (7,7') gener- 
ated by reshuffling intra-encounter connections of loops 
within an Z-encounter with I > 2 or within several en- 
counters; see Fig. lb. In general, such reconnections arise 
in vi Z-encounters, I = 2, 3, . . .; we collect these multiplic- 
ities in a vector v. The total number of encounters is 
V = X)i=2 v i wn il e L — X)z=2 gi yes t ne number of the 
encounter stretches; the total number of loops is L + 1. 
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Each loop is represented by one loop-time integral as a 
factor and each encounter by an (I — l)fold (s, u)-integral. 

N ,n 

The (s, w)-integrals are done as above, e T « cnc — > 



T^Cc and (^Jdsdue 1 



su/h\ 



L-V 



T. 



V-L 



H 



In the last member of (0 the sum over the possible 
vectors v has been split into two, one with v restricted 
as L — V = n and a subsequent one over the integer n. 
The partial sum with fixed n can be carried out with the 
help of a recursion relation for the Af(v) obtained in Q 
and reads (see the first technical note below) 



(£-V=n) 

£ (-l) v N(v) = (l-2/f3Y 



(8) 



Upon inserting that partial sum into (JJJ we recover the 
full l/7V-expansion (JTJ of the conductance, now for a sin- 
gle device rather from an ensemble average. 

The transition between the orthogonal and unitary 
symmetry class in a weak magnetic field has been dis- 
cussed in Ref. 0, 0, 0- Saito and Nagao [l3] have 
recently accounted for a magnetic field in the semiclas- 
sical r-expansion of the spectral form factor, up to r 3 . 
We here extend their ideas to the 1/iV-expansion of 
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the conductance. The additional action accumulated in 
the magnetic field B on some trajectory interval CD is 
8cd = f Jq Avdt with A the vector potential. If the 
dynamics is chaotic the effect of the field has been shown 
to be representable by Gaussian white noise Q which 



yields the average (e %26cD l h ') Q = e btCD - here b, scaling 

like 9c D , is (B/h) 2 times a system specific parameter. 
We shall see presently that net contributions of that type 
come both from loops and encounter stretches. 

We consider a pair (7, 7') and denote by 6i the mag- 
netic action gained on the z-th loop with duration tj, in 
7. If that loop is passed in the same direction by 7 and 
7' that action cancels in (S 1 — S 7 ')/h; otherwise it is 
doubled. Defining a quantity Si equal to and 1 in 
the first and second case, respectively, we can write as 
(e l2<5i8i / ri ) G = e~ b5iti the additional factor in the contri- 
bution of the pair to the conductance associated with the 
loop H" . Now consider an encounter "a" and denote by 
9 a the magnetic action accumulated on a single stretch of 
this encounter passed in some ("positive") direction. The 
magnetic action gained on all stretches of "a" in 7 will be 
a a 9 a with <r a the difference between the number of the 
encounter stretches traversed in the positive and nega- 
tive direction. The corresponding action for the partner 
7' will be <j' a a ; generally a' a ^ a a . Averaging yields 
the factor (Q l2 ^a6 a /hy _ e -iJ, a bt enc j n ^ e contribution 
of the orbit pair, with the integer /i a = (a a — a' a )/2. 

We shall parametrize the magnetic field by £ = 
bT H /N oc (B/h) 2 and denote by j/ = Y%=i S i the 
number of loops whose sense of traversal is changed in 
7' relative to 7. Finally, we refine our definition of 
the number of structures by introducing the numbers 
J\f(v, v, {fJ-a}) of topologically different trajectory pairs 
with fixed v, v, {/i Q }. The above result Q is then refor- 
mulated for the magnetoconductance, 

" {^-Ml) (9) 
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N 
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v,v,{n a } 



£ + 1 i* 00 

M(v, u, {^a})TT / dUe 

1 Jo 



-(l+£<5«)ti. 




n=l v,u,{^i a } 



(-i)^>,K})nLi(i+^) 



Obviously, the limiting case £ — > yields the orthogonal- 
class conductance Q due to Ylu{^} "* N(v, v, {^ a }) — 
J\f(v); the opposite limit £ — > 00 is less transparent. The 
coefficients of the resulting 1/iV-cxpansion for low or- 
ders n are easily obtained by checking the contributing 
structures of orbit pairs for their values of v, {fi a } and 
counting the multiplicities Af(v, we get 



Tr(^) 
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(10) 



1 £ 4 + 4£ 3 + 13£ 2 + 2£+l 



(i + E 



It is seen that when £ — > 00 we do recover the unitary- 
class conductance. Order by order, our expression ltTU|) 
coincides with the RMT-based expansion of the integral 
obtained in [9j. For n > 3, it appears that our semiclas- 
sical sum © is easier to evaluate than the corresponding 
term in . The expansion coefficients up to order n = 8 
are listed in the second technical note below. 

Technical note I: The pairs of open trajectories of 
interest here and pairs of periodic orbits with the same v 
have the same number Af(v) of structures. This is a mat- 
ter of definition: we had associated in Q exactly as many 
structures with a single periodic orbit pair, as there are 
non-equivalent trajectory pairs produced by cutting one 
of the orbit loops. We may thus use the recurrence rela- 
tions for the number of periodic orbit structures J\f(v) to 
find the sum JSJ • We need the special case incorporating 
equations (44,60,30) of the online version of 0, 
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Niy) 



Z 7^3 T E( fc + 1 K + i^') 

V ' k>2 



(J 



N{v") 



(11) 



Here v' (denoted y^ k ' 2 ^ k+1 ^ in Q) is the vector obtained 
from v by decreasing the components vi and i>& by 1 
and increasing Vk+i by 1; the vector v" (v^ 2 ^ of Q) 
is obtained by decreasing V2 by 1. We multiply both 
parts of by (— l)^) and sum over all v with fixed 
difference L(v) — V(v) — n. In the double sum over v, k 
we replace v by v ' as the summation variable and use 
L(v') = L(v) - l,V(v') = V(v) - 1; we then drop the 
prime in v'. Similarly, we switch in the right hand side 
to v" as the summation variable, taking into account 
L(v") = L(v) - 2, V(v") = V(v) - 1, L(v") - V(v") = 
n — 1. Dropping the double prime we obtain 



L{v)-V{v)=n 

E 



m 



2v 2 + 2~^{k + l)v k+1 



k>2 



(-i) y ^(«). 



(12) 



The expression in the square brackets is by definition 
L(v) and thus cancels with the denominator. The result- 
ing recursion relation S n = (1 — 2/ f3)S n —i for the sum 

L-V=n 

S n = E (-l) v M(v) is solved by S n = S (l - 2/(3) n . 

v 

We conclude S n = for all n > 1 in the unitary case. For 
the orthogonal case we have to fix the constant factor So. 
Now for n = 1 we have but one 2-encounter (d 2 = 1, all 
other components of v zero) and a single structure, i.e. 
J\f(v) — 1 such that Si = — 1 and thus So = 1. 
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Technical note II: We here extend the 1/N expan- 
sion 11U|) of the conductance in the presence of a magnetic 
field. Defining the expansion coefficients as 

(13) 

n=Q 

and introducing the abbreviations <7fe(£) = 1 + 
fc 2 £, p(£) = 1/(1 + £) = q^ 1 we have, by numerically 
determining the numbers Af(v, v, {/^a}), 

C = 1, C 1 = -p, C 2 = P 2 , 

C 3 = -p + qsP 4 -qlp\ 

C 4 = p 2 -p i ~4q 3 p 5 + 5q 2 p 6 7 

C 5 = -p + q 3 p 4 - q%P 5 + 8q 5 p 6 - 24q 2 q 4 p 7 

-I2q 2 p 7 + 49q 2 q 3 p s -2lq 4 2 p 9 , 
C 6 = p 2 - p i - 4g 3 p 5 + 5(? 2 / + 4p 6 + 16qsp 7 

-A8q 5 p 7 + I68q 2 q4p 8 - 2\q\p 8 + 84gfp 8 

-392q 2 q 3 p 9 + 189<7 2 y°, 
C 7 = -P + A 3 -P 5 ?2+ 8/(75 - 12/g 3 2 - 24A294 

+I8OA7 + 49p 8 9 2 g 3 - 21p 9 9 | - 608p 9 g3 g5 



-276AI - 720A 2 (7 6 + 1728p 10 g 2 g5 
+2952p 1( W 3g4 + 464p 10 g 3 - 3240p n ^ g4 
-4440p n (7 2 2 (7 3 2 + 5445p 12 (7 2 4 (7 3 - 1485p 13 (7 2 6 , 
C 8 = p 2 +4p 6 -p 4 -36p s + 2l6p w q 2 ~21p 8 q 2 
-1485p 12 ^ + 5A2 + ^p w qt 
+ 19305p 14 g 2 3 - 160pV 3 + lQp 7 q 3 
+2880p n q 2 2 q 3 - 4p 5 q 3 - 392p 9 q 2 2 q 3 
-65340p 13 g 2 g 3 - 1080p w q 2 q4 + 168p s q 2 q 4 
+35640p 12 g| 94 - 624p 1Q q 2 + 84p 8 ql 
+48840p 12 g|g 2 + 288p 9 c?5 - 48p 7 95 
-17280p n g| g5 - 29520p n q 2 q 3 q 4 - 4640p n 9 | 
+6480p w q 2 q 6 + 2484p w q 2 + 5472p w q 3 q 5 
-1440pV- 

A recursive treatment of all orders in the expansion 113|) 
is called for. 

We are indebted to Dmitry Savin for useful discus- 
sions and the Sonderforschungsbereich SFB/TR12 of the 
Deutsche Forschungsgemeinschaft for financial support. 
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